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I study the properties of the scalar a-meson [also referred to as /o(600)] at nonzero temperature 
. in the 0(N)-mode\ in the framework of the Cornwall- Jackiw-Tomboulis formalism. In the standard 

■ Hartree (or large-N) approximation one only takes into account the double-bubble diagrams in the 

Pj^ ' effective potential. I improve this approximation by additionally taking into account the sunset 

diagrams, which leads to 4-momentum dependent real and imaginary parts of the Dyson- Schwinger 
equations. By solving these and the equation for the chiral condensate selfconsistently, one gets the 
decay width and the spectral density of the cr-meson. I compare the results in the case with explicit 
chiral symmetry breaking with the chiral limit. I found, that the 4-momentum dependent real part 
of the selfenergy do not lead to qualitative changes in the spectral density. 
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In the low-energy regime of QCD, chiral symmetry is spontaneously broken by a nonvanishing 
quark condensate (qq) ~ (300MeV) 3 . (In nature this symmetry is also explicitly broken by a nonvan- 
ishing quark mass.) At high temperatures, T ~ (qq), one expects a "melting" of the condensate and 
therefore a (partial) restoration of this symmetry. A major goal in modern hadron physics is to find 
clear indications for this chiral restoration. A reasonably way is to study the properties of strongly 
interacting matter near the critical temperature of this transition, e.g. in heavy ion collisions. 

In this work I will focus attention on the sector of light mesonic degrees of freedom at finite 
temperature. An adequate field-theoretical framework for this is the linear sigma model with O(N) 
^ . symmetry, which has been extensively studied in the last decades, cf. P, IE IE IE IE IE IE IE EE 

fill \13L Il4j and citations therein. In this model the fundamental degrees of freedom of QCD (the 
quarks and gluons) are integrated out, and new (effective) degrees of freedom (the mesons) arise, 
which in the O(N) model are the scalar cr-meson and the (N — 1) pseudoscalar pions. These are 
chiral partners and thus become degenerate in mass in the chirally restored phase. (For an extension 
to finite densities, cf. e.g. [lEIlEEl) 

In general, neither QCD nor the O(N) model can be solved analytically. A standard way to get 
results is ordinary perturbation theory which is an expansion up to a certain order in the coupling 
constant g. Unfortunately, this scheme breaks down at nonzero temperature [l^. The reason is that 
the new energy scale introduced by the temperature can conspire with the typical momentum scale 
p of a process so that gT/p is no longer of order g, but can be of order 1 |l9|,|20j. For this reason, 
one has to find a proper way to resum whole classes of diagrams of order gT/p. 

The resummation method I use here is the so called Cornwall- Jackiw-Tomboulis (CJT) formalism 
|2l| at finite temperature. This scheme is equivalent to the ^-functional approach of Luttinger and 
Ward [2E] and Baym (2^. In the CJT formalism (assuming translation invariance), the Dyson- 
Schwinger and the condensate equations are given by the minimum of an effective potential 

V[a,G] = U(a) + \ f baG-^ + i f [G'^Q, a)G(Q) - 1] + V 2 [a, G] (1) 

1 JQ 1 JQ 

where cf and G are the expectation value of the one and two-point function in the presence of external 
sources, U(a) is the tree-level potential, G _1 the inverse tree-level two-point function, and V<2,[cr, G] 
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the sum of all two-particle irreducible (2PI) vacuum diagrams with internal lines given by G (the 
definition of Jq ■ ■ ■ is given at the end of this introduction). Taking into account all 2PI diagrams 
in the functional V2 would be equivalent to solve the exact theory, which is impossible. To get 
numerical results one has to truncate this sum. 

In the standard Hartree (or Hartree-Fock) approximation of the O(N) model, one only takes 
into account the double-bubble diagrams, shown in Figs. ^ a i b, c. In these approximations the 
selfenergies of the particles are only real valued, therefore no nonzero decay width effects are included. 
The difference between the Hartree and the Hartree-Fock approximations is that in the Hartree (or 
large- N) approximation all terms of order ~ 1/N are neglected on the level of the Dyson-Schwinger 
and the condensate equations. 

In |24| we presented the so-called improved Hartree-Fock approximation, which takes additionally 
into account sunset type diagrams, shown in Figs. Ed, e. This leads to 4-momentum dependent 
real and imaginary parts in the Dyson-Schwinger equations. Indeed, in that paper, we neglect 
the 4-momentum dependent real parts of the Dyson-Schwinger equations for simplicity. In the 
present work, I include them, although just in the Hartree approximation. This leads to a vanishing 
imaginary part of the pion selfenergy (because it is of order ~ 1/N), i.e., the spectral density of the 
pion is a delta function. For the cr-meson the contribution from the a — » 2a decay is also of order 
~ 1/N, and vanishes, but the contribution from the a — > 27r decay remains. Therefore, the spectral 
density has a nonzero width, as expected for the cr-meson with a very large vacuum decay width, of 
T a = (600 — 1000) MeV [2q. In the following, I call this the improved Hartree approximation. 

The paper is organized as follows. In Sec. II I derive the Dyson-Schwinger and the condensate 
equations of the O(N) model in the improved Hartree approximation. In Sec. Ill the numerical 
results for the N = 4 case are presented in the chiral limit and in the case with explicit chiral 
symmetry breaking for the parameter sets summarized in Tab. [fl Section IV concludes the paper 
with a summary and an outlook. The technical derivation of the integrals, which appear in the 
Dyson-Schwinger and the condensate equations are given in the appendix. 

I denote 4- vectors by capital letters, X = (xo,x), with x being a 3- vector of modulus x = |x|. 
I use the imaginary time formalism to compute quantities at finite temperature. Integrals over 
4-momentum K = (feo,k) are denoted as 



where T is the temperature and uj n = 2nnT, n = 0, ±1, ±2, . . . are the bosonic Matsubara frequen- 
cies. Units are % = c = ks = 1- The metric tensor is = diag(+, — , — , — ). 

II. THE O(N) MODEL IN THE IMPROVED HARTREE APPROXIMATION 

In this section I discuss the including of the sunset diagrams in the standard Hartree approxima- 
tion. The Lagrangian of the O(N) linear sigma model is given by 



where <f> = (0i,7r), with the first component <j>\ corresponding to the scalar cr-meson [also referred 
to as /o(600)] and the other components tt — {<j>2, 4>n) corresponding to the pseudoscalar pions. 
For H = and /i 2 > 0, the Lagrangian is invariant under O(N) rotations of the fields. For negative 
values of the bare mass /i 2 , the symmetry is spontaneously broken from 0(N) to 0(N — 1), which 
leads to N— 1 Goldstone bosons (the pions). The parameter H breaks the O(N) symmetry explicitly 
and give a mass to the pions. 

The parameter H is given as a function of the vacuum mass m n , and the vacuum decay constant 
f„ of the pion : H = m^f^. In this work I compare the chiral limit, m, = and / w = 90 MeV, 




(2) 



C(4>) = - (d,<j> ■ d»4> - /i 2 • 0) - j^(4> ■ 4>) 2 + Hfa 



(3) 
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cr-meson 
vacuum mass 


explicit chiral symmetry breaking 




chiral limit 


m n = 139.5 MeV, f„ = 92.4 MeV 




= OMeV, /, = 90 MeV 


H = (121.6 MeV) 3 




H = 


m a = 400 MeV 


A = 8.230 


A = 


9.877 




fjj 2 = -(225.41 MeV) 2 


M 2 = 


= -(282.84 MeV) 2 


m a = 600 MeV 


A = 19.043 


A = 


22.222 




/i 2 = -(388.34 MeV) 2 


M 2 = 


= -(424.264 MeV) 2 


m a = 800 MeV 


A = 36.341 


A = 


39.506 




fj 2 = -(539.27MeV) 2 


M 2 = 


= -(565.685 MeV) 2 



TABLE I: The masses and decay constants at vanishing temperature and the corresponding parameter sets 
for the 0(4) linear sigma model for the two symmetry breaking patterns studied here. 



with the case of explicit chiral symmetry breaking, — 139.5 MeV and = 92.4 MeV. The 
bare mass /i 2 and the four-point coupling A depend additionally on the vacuum mass of the cr- 
meson 77i CT : [x 2 = — {m 2 — 3m 2 )/2, and A = N(m 2 — m 2 )/(8/ 2 ). The decay width of the cr-meson 
in vacuum is very large, Y a = (600 — 1000) MeV, therefore the vacuum mass is not well defined, 
to ct = (400-1200) MeV 25]. I compare the results for m a = 400, 600, and 800 MeV. The parameter 
sets of the model, for N = 4, are summarized in table U 

The effective potential for the 0(N) model in the CJT formalism, Eq. 0J, is given by [fl fill] 

V[a,S,P] = l^ + ^-Ha + i f [hiS'^Q) + S~\Q; a) S(Q) - 1 ] 

■/ [In P- 1 (Q) + P-\Q; a) P(Q)-1] +V 2 [d,S,P] , (4) 
JQ 



N — 1 



where cr, S, P are the expectation values of the one- and two-point functions in the presence of 
external sources (I omit the dependency on the pseudoscalar one-point function n, because the 
vacuum has even parity, and therefore ff = 0). At tree-level the expectation value of a can be 
calculated analytically 



2 e 

<7 = u = ^-^Tv^ C0S ^ 



HN ( 12A x 3/2 



8A V A^/i 2 



(5) 



The quantities S 1 and P 1 are the inverse tree-level propagators for scalar and pseudoscalar 
mesons, 

S-\K:a) = -K 2 +ml(a) , p-\K- a) = -K 2 + m 2 (a) , (6) 
where the tree-level masses are 

2/-\ 2 i 12 A _ 2 2/-\ 2 i ^A -2 m\ 

As discussed in the introduction, in the functional V2, additionally to the three double-bubble 
diagrams shown in Figs. ^ a j b, c, the two sunset diagrams shown in Figs.^d, e, are taken into 
account, 
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FIG. 1: The set of two-particle irreducible diagrams included in the effective potential of the CJT-formalism. 
The diagrams a,b, and c are the double-bubble diagrams, and d and e are the sunset diagrams. Full lines 
are scalar particles, dashed lines are pseudoscalar particles. The four-particle vertex ~ A is represented by 
a full squared and the three-particle vertex ~ Act is represented by a full triangle. 



V 2 [a,S,P] = (JV+1)(JV-1) — 



+ 2(N-1) 



P(Q) 
§(Q) I P(Q) 
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S(Q) 



\ N J 



S(L)P(Q)P(L + Q) 



-3!( — ) / / S(L)S(Q)S(L + Q). 



(8) 



To get the expectation values for the one- and two-point functions in the absence of external sources, 
a, S, and V, one has to find the stationary points of the effective potential Minimization of the 
effective potential with respect to the one-point function SV/Sa = 0, leads to an equation for the 
scalar condensate er, 

4A 4A f 

H = 3 + A / [3S{Q) + (N _ 1)Vm 

a 2(N-1) [ [ S(L)V(Q)V(L + Q) +3! f [ S(L)S(Q)S(L + Q) . (9) 
J l Jo Jl Jo 



Using the fact that a 2 ~ N [cf. Eq. JSJ], this equation becomes in the large- N limit, 

H = aL 2 + ^ a 2 +N j V(Q) 



(10) 



The minimization of the effective potential with respect to the two-point functions, SV/5S — and 
SV/SP = 0, leads to the Dyson-Schwinger equations for the scalar and pseudoscalar propagators, S 
and Q, 



S-^K;*) = 5 , " 1 (A';cr) + S(if;ff) , V~ X {K-. a) = P~ X {K\ a) + U(K; a) 
Here I introduced the selfenergy of the scalar 
4A 



(11) 



N 



/4Acr\ 



Q 
2 r 



S(Q) +(N-1) / V(Q) 

JQ 

2(N -I) I V(K - Q)V(Q) +3-3! f S(K - Q)S(Q) 
Jo Jo 



(12a) 
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and pseudoscalar fields 



U(K:a) _ 



J S(Q) + (N + 1) J V(Q) + (jf^j 4 J V(K - Q)S(Q) . (12b) 

In the large- N limit, all cr-meson tadpoles and a part of the pion tadpole (~ J S and ~ f V) vanish, 
and only a part of the pion term (~ J VP) remains, 

V{K-a) = ^N J V(Q) + 2Nj V(K-Q)V(Q), IL(K;a) = ^N J V{Q). (13) 



IQ \ JV / JQ 

The tadpole contributions have no imaginary parts, therefore 
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JmE(K;a) = (^jf) 2Nlm J V(K - Q)V(Q) , Imll 



(14) 



The whole imaginary part of the cr-meson selfenergy depends on the 4- momentum vector, K — (k, u>), 
but the real part can be split into a term which do not depend on K, [ReE]i, stemming from the 
tadpole contribution, and a term which is 4-momentum dependent, [Re S(if, f)]2, stemming from 
the sunset diagram, 

ReZ(K;a) = [Re E] i + [Re E(JT; cr)] 2 , (15a) 

where 

[ReE] x = ^iV J v(Q), [Rel2(K;a)} 2 = 2 NReJ^V(K- Q)V{Q). (15b) 

As shown above, the 4-momentum dependent terms in the pion selfenergy vanish in the large-iV 
limit, and the 4-momentum independent term is the same as for the cr-meson, Re II = [ReE]i. 

In this work I want to calculate the spectral density of the cr-meson p a . To this aim, I rewrite the 
Dyson-Schwinger equations (flip and the equation for the chiral condensate I|1Q(I as functions of the 
spectral density of the pion p^. In the large- N limit, the imaginary part of the pion vanishes, and 
therefore the spectral density is just a delta function 

7T 

p„(w, k) = — — {5[lo - io-n{k)} - 5[ui + oj^(k)]} , (16) 



with support on the quasiparticle energy for the pion u w (k) = \/k 2 + M 2 (cr), where I defined a 
4-momentum independent mass for the pion M n , 

Mj(cr) = ml(<r) +ReU . (17) 

In the chirally broken phase (cr 7^ 0) the imaginary part of the cr-meson is nonzero, therefore the 
spectral density assumes the following form: 

/ .\ 2ImE(u;,k; a) 

Pa{UJ ' ' [^ 2 -fc 2 -TO2(CT)-ReE(w,k;CT)] 2 + [ImE(w,k;cr)] 2 ' 1 ' 

where the quasiparticle energy for the cr-meson w CT (k), is given by the solution of u;^(k) — k 2 — 
m 2 (cr) — Re E[w (T (k), k; a] = 0. In the chirally restored phase (cr = 0) the 4-momentum dependent 
real and imaginary parts of the a-meson selfenergy vanish, therefore also the spectral density of the 
cr-meson becomes a delta function, 

p a (u>, k) = {S[u> - w a {k)] - 5[uj + u a {k)]} , (19) 

uj a {k) 



() 



with support on the quasiparticle energy for the cr-meson u> a (k) = \fk 2 + M% (a) , where I defined a 
4- momentum independent mass for the cr-meson M a , 

M^a) = ml(a) + [ReT,} 1 . (20) 

The selfenergies can be written as functions of the pion spectral density, exclusively, which is just 
a delta function as discussed above. Note, that for the real parts, I only consider temperature- 
dependent contributions and neglect the (ultraviolet divergent) vacuum parts, which is a simple way 
to renormalize the integrals (the imaginary parts are not ultraviolet divergent). The 4- momentum 
independent term is just the standard tadpole integral [Tlj . 



1 f°° 

Ren = [Rcs]x = — / dqc/M^rVM*?)]. 

2tt 2 Jo 



(21) 



The derivation of the equations for the 4-momentum dependent terms is given in the appendix, 

— J g-jj / dq qlu^q^Q (\q - q\ < k < q + q) 

x {1 + /K(gd)] + /K («)]}, (22a) 

— J ^2^ P / d< 7l<7l d< 72<72 6 (kl - g2| < fc < 91 +<?2) [^ 7 r('7l)w7 r ('?2)] _1 

/K(gi)] + /K(<&)] , /K(gi)] + /K(g 2 )] 



-/K(gi)] + /K(g 2 )] | -/K(gi)] + /K(g 2 )] 
w w (<jfi) - 0J„(q2) - co uj^iqi) - u 7r (q 2 ) + uj 



(22b) 



where /(w) = l/[exp(w/T) — 1] is the Bosc-Einstein distribution function, go = vP^ — ^(g)] 2 — Af^, 
and P J . . . denotes the principal value of the integral. Note that O (|gi — e/ 2 | < fc < gi + g 2 ) /fc = 
2<5(c/i — g 2 ) in the limit k — > 0, which can be used to perform the g-integration, cf. Eqs. (|33|) . and 
(I39|l . An appropriate way to perform the principal value numerically is discussed in the appendix, 
cf. Eqs. and 

The spectral densities have to obey a sum rule |2(J, J*^° dwo;p a - i7r (a;,k) = 2n , which is a priori 
fulfilled for the pion spectral density (|16fl . because of the normalization. As mentioned above, 
this is not the case for the cr-meson in the chirally broken phase. The main reason for a possible 
violation of the sum rule is due to neglecting terms of the order ~ 1/N in the selfenergy, which 
leads to a loss of spectral strength. (I found that the inclusion of the 4-momentum dependent real 
part [ReS] 2 , is close to negligible for the validity of the sum rule.) Other possible problems arise 
from the numerical realisation of the spectral density on a finite energy-momentum grid, which I 
minimized by using a large and fine grid. As discussed in [24[, if the sum rule is not fulfilled, I use 
the following way to restore it. If the imaginary part of the cr-meson selfenergy is very small (i.e. 
the spectral density is narrow), I add a numerical realisation of a delta function S nnm to the spectral 
density: p' a (ui, k) — > po-(w, k) + c\ ■ <5 nU m[^ — w CT (fe)]. On the other hand, if the imaginary part of the 
selfenergy is large enough (compared with the lattice spacing) I just multiply the spectral density by 
a factor: p' a (uj, k) — > c 2 ■ Po-(to, k). The constants c\ and c 2 are adjusted in that way, that p' a fullfills 
the sum-rule on the (finite) energy-momentum grid. (I checked, by comparing the results with the 
case c 2 = 1, that this procedure does not lead to major quantitative changes.) 

The numerically scheme for the improved Hartree approximation is the following. At first, 
one solves the standard Hartree approximation, i.e., the condensate and the Dyson-Schwinger 
equations, Eqs. (|10f> and ^ X 1 p . with [ReS]2 = [ImS]2 = 0, to get the chiral condensate and 
the 4-momentum independent mass of the pion, a and M n . With these results, one calcu- 
lates the 4-momentum dependent real and imaginary parts Q22H of the cr-meson selfenergy, and 
finally the spectral density, Eq. I|18fl or CHJ- Then, the decay width of the cr-meson T a is 
defined as [2(| [2?]]: E (T (A;) = lmT,[uj cr (k),k;a]/uj cr (k), where uj a {k) is given by the solution of 
colik) - k 2 - ml(a) - ReE[w CT (fc), k; a] =0. 
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FIG. 2: The chiral condensate a (left column) and the 4-momentum independent mass of the cr-meson and 
the pion, M a and M w (right column), as function of temperature T and cr-meson vacuum mass m a . In the 
upper row are shown the results with explicit chiral symmetry breaking and in the lower row the results in 
the chiral limit. 



III. RESULTS 



In this section I present the numerical results for the linear sigma model with 0(4) symmetry in 
the improved Hartree approximation, for the parameter sets given in table IJ 



A. The mass 



I start the discussion of the results with the chiral condensate a and the 4-momentum independent 
mass of the cr-meson and the pion, M a and M ff , as defined in Eqs. (|20|l and l|17|l . Note, that the 
4-momentum independent mass M a is generally not the physical (pole) mass of the cr-meson, because 
the pole of the cr-meson spectral density is additionally shifted by the 4-momentum dependent part 
of the selfenergy [Re Y,(K : cr)] 2 , cf. Fig.0 

In the upper row of Fig.[5]the results for the chiral condensate a (left column) and the 4-momentum 
independent mass for the pion and the cr-meson, M T and M a (right column), are shown as functions 
of T and m a , in the case of explicit chiral symmetry breaking. The results show the behaviour of 
a crossover transition. Neither the chiral condensate nor the 4-momentum independent mass of the 
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FIG. 3: The 4-momentum dependent real part of the cr-meson selfenergy [ReE]2, at a fixed momentum of 
k = 1 GeV, for the case with explicit chiral symmetry breaking (full line) and the corresponding result for 
the chiral limit (dotted line) at a temperature of T = 120 MeV. 



scalar particle M a become equal zero, even for high temperatures. Thus, M n and M CT , become only 
approximatively degenerate for high temperatures. In the lower row of Fig. [21 the corresponding 
results for the chiral limit are presented. The results show the behaviour of a second-order phase 
transition, which agrees with the predictions made by Pisarski and Wilczek [28|. The condensate a 
becomes equal zero at a critical temperature T x , therefore M a = M w for T > T x . The condensate 
(accordingly T x ) does not depend on the vacuum mass m a , which can be understood as a consequence 
of the condensate equation l|10|l in the chiral limit (H = 0), 



= + 



4A 
N 



N I P(Q) 
IQ 



For = the integral can be performed analytically 



2 



4A 
— c 

N 



12 



T 2 N 
12 



(23) 



(24) 



[cf. Eq. 0], and the critical temperature is given by T x = ^jYljNj-,, w 160 MeV for the 0(4) 
model, which agrees with the numerical results. 

The 4-momentum dependent real part of the a- meson selfenergy [ReS]2 is shown in Fig. [21 at 
fixed momentum of k = 1 GeV as a function of energy. It is larger in the chiral limit as in the 
case with explicit chiral symmetry breaking, because [Re£] 2 ~ A 2 ~ (m 2 — to 2 ) 2 is maximal for 
to 2 = 0. Note that [ReS]2 is small compared to the (squared) 4-momentum independent mass of 
the cr-meson, shown in Fig. [21 



B. The decay width 

After calculating the (4-momentum dependent) imaginary part of the cr-meson selfenergy (|14|) . the 
decay width T a is calculated at the quasiparticle energy: T a (k) — Im E[w CT (fc), k; a\/ui a (k), where 
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FIG. 4: The decay width of the a- meson r CT , as a function of temperature T and a- meson vacuum mass m a , 
at momentum k = 1 GeV, in the case with explicit chiral symmetry breaking (left) and in the chiral limit 
(right). 



uj a (k) is given by the solution of w 2 (k) — k 2 — m 2 (a) — Re Y\u) a (k), k; a] = 0. The imaginary part of 
the sclfcnergy, and therefore the decay width, of the pion is of order ~ 1/N, and therefore neglected. 

In Fig. 0] the decay width of the cr-meson T a is shown as a function of temperature T. The 
qualitative behaviour is similar in all cases, but the decay width is larger in the chiral limit compared 
to the case with explicit chiral symmetry breaking, because ImE ~ A 2 ~ (to 2 — m 2 ) 2 is maximal 
for to 2 — 0, cf. above discussion of [ReS]2- The decay width is a strictly monotonic decreasing 
function with temperature, and becomes approximatively zero (equal zero) for temperatures larger 
than ~ 200 MeV (T x ~ 160 MeV) in the case of explicit chiral symmetry breaking (the chiral limit). 
This is a consequence of the (partial) restoration of the chiral symmetry, the masses of the chiral 
partners become (approximatively) degenerate in mass, and therefore the phase space of the a — ► 2ir 
decay is squeezed. The dependence on the vacuum mass of the cr-meson is significant. The reason 
for this is that T a ~ ImE ~ m^., which agrees reasonably with the results, ImS = 1 : 3 : 16 for 
to ct = 400 : 600 : 800 MeV. 

In Fig. the decay width of the cr-meson in the improved Hartree (triangles) is compared with 
the improved Hartree-Fock [24| (squares) approximation. The main difference stems from the com- 
binatorial factor in front of the two-pion term (~ / VV) in the imaginary part of the selfenergy [cf. 
Eqs. H12af) with l|14l) ]. A part of this contribution is of order ~ 1/N, and neglected in the Hartree 
approximation. Therefore, this factor is 2 ■ N — 8 in the improved Hartree, and 2 ■ (N — 1) = 6 in 
the improved Hartree-Fock approximation, which would lead to a decay width of an factor w 1.33 
larger in the improved Hartree approximation. The remaining difference, shown in the plot, comes 
from the cr-meson term (~ 3-3! J SS) in Eq. I|14|) . which vanishes in the large- N limit. 



C. The spectral density 



Finally, after solving the coupled condensate and Dyson-Schwinger equations selfconsistently, 
Eqs. HI U| l and i|ll|) . the spectral density of the cr-meson p a is given in the chiral broken phase 
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FIG. 5: The decay width of the cr-meson F a , in the improved Hartree (triangles) and in the improved 
Hartree-Fock (squares) approximation, as a function of temperature T, at momentum k — 1005 MeV and 
cr-meson vacuum mass of m a = 600 MeV, in the case with explicit chiral symmetry breaking. 



[a ^ 0) by Eq. iJTgfl. and in the restored phase (a = 0) by Eq. (|T§1) . 

In Fig. the spectral density is compared for all possible parameter sets (cf. table I) at a 
fixed momentum of k = 1 GeV. The spectral density does not exhibit a pronounced peak at the 
quasiparticle energy because this energy is large enough for the decay into two pions. As 

discussed above, the decay width of this process is large and becomes even larger for increasing m CT , 
as shown in the plot. A remarkable difference between the chiral limit (right) and the case with 
explicit chiral symmetry breaking (left) is, that p a ~ around uj « 1 GeV in the case with explicit 
symmetry breaking but not in the chiral limit. The reason is, that the threshold energy for cr — ► 2n 
is lo = 2M W , and therefore zero in the chiral limit. For lu < k the cr-meson is Landau-damped. 

To illustrate the effects of neglecting terms of order ~ 1 /N in the improved Hartree approximation, 
I compare in Fig. {7\ the results with the improved Hartree-Fock [24J approximation, in the low- 
temperature regime at T — 160 MeV (left), and in the high-temperature regime T = 320 MeV 
(right). As show, the terms of order ~ 1/N (included in the improved Hartree-Fock approximation) 
leads to a nonzero decay width of the pion, and therefore to a washed-out cr-meson spectral density. 
One aim of this work was to study the influence of the 4-momentum dependent real parts of the 
selfenergy. I found, that they does not lead to qualitative changes of the spectral density. To get an 
quantitative estimate, I consider the relative change between the spectral density with and without 
this real part, averaged over the energy-momentum grid, which is rather small, e.g. for the results 
shown in Fig. □ 5.61 ± 2.90% for T = 160 MeV, and 0.44 ± 0.01% for T = 320 MeV. 



IV. CONCLUSIONS 



In this paper I systematically improved the standard Hartree (or large- N) approximation of the 
O(N) linear sigma model by taking into account, additionally to the double-bubble diagrams (shown 
in Figs, n a, b, c), the sunset diagrams (shown in Figs. Ed, e) in the 2PI effective potential of 
the CJT formalism. This leads to 4-momentum dependent real and imaginary parts of the Dyson- 
Schwinger equations. I solve these and the equation of the condensate selfconsistcntly to get the 
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FIG. 6: The spectral density of the <r-meson p CT , as a function of energy u), and cr-meson vacuum mass m a , 
at a momentum of k — 1 GeV and a temperature of T = 100 MeV, in the case with explicit chiral symmetry 
breaking (left) and the chiral limit (right). 

decay width and the spectral density of the cr-meson, for the parameter sets summarized in table I. 

First, I presented the results for the real parts of the Dyson-Schwinger, and the condensate 
equations. The 4-momentum independent parts, i.e., the 4-momentum independent masses and the 
chiral condensate exhibit a crossover transition in the case with explicit chiral symmetry breaking 
and a second order phase transition in the chiral limit, as predicted by Pisarski and Wilczek 28] . The 
4-momentum dependent real parts are small compared to the (squared) 4-momcntum independent 
masses. 

The decay width, presented in the second part of the results, show qualitatively the same be- 
haviour in all cases. It is a strictly monotonic decreasing function with temperature, and becomes 
approximatively zero (equal zero) for temperatures larger than ~ 200 MeV (T x w 160 MeV) in the 
case of explicit chiral symmetry breaking (the chiral limit) . This is a consequence of the (partial) 
restoration of the chiral symmetry, the masses of the chiral partners becomes (approximatively) 
degenerate in mass, and therefore the phase space of the a — > 2ir decay is squeezed. The dependence 
on the vacuum mass of the cr-meson is significant, because of r CT ~ ImE ~ m\. 

Finally, the behaviour of the decay width and the real parts are reflected in the spectral density. 
In the low-temperature regime, the spectral density is a very broad function in energy (and becomes 
even broader for larger vacuum mass m a ) due to the cr — > 2 7r decay and becomes more and more a 
delta function for increasing temperatures. The 4-momentum dependent real part of the selfenergy 
do not lead to qualitative changes. 

The next step is to include also terms of order ~ 1/N, which would lead to a nonvanishing decay 
width of the pion and therefore to washed-out spectral densities. Other possible future pro jects are 
the inclusion of more degrees of freedom, for instance baryons [2t|, and vector mesons |30i l3ll l32| . 
The latter are of particular importance, since in-medium changes in the spectral properties of vector 
mesons are reflected in the dilepton spectrum |31| which, in turn, is experimentally observable in 
heavy-ion collisions at GSI-SIS, CERN-SPS and BNL-RHIC energies. 
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FIG. 7: The spectral density of the cr-meson p a , as a function of energy u>, at a momentum of k — 1 GeV, and 
a vacuum cr-meson mass of m a = 600MeV, in the case with explicit chiral symmetry breaking. The results 
are shown for the improved Hartree (full line) and the improved Hartree-Fock (dotted line) approximation 
at temperatures of T = 160 MeV (left), and 320 MeV (right). 
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Appendix: The 4-momentum dependent parts of the selfenergy 



In this section I derive the equations for the 4-momentum dependent imaginary and real part 
of the selfenergy (stemming from the sunset diagram). As discussed in |2l| this diagram can be 
expressed as a function of the spectral density p v in the imaginary time formalism 



E(-iuj m ,k) 



2N 1 



\~N~) (2tt) 4 k J_ 



du>i duj\ 



1 + /M + /M 

iuj rn + UJ\ + L02 



dqi qi dq 2 q 2 6 (\qi - q 2 \ < k < q 1 + q 2 ) Pn(^l,Ql)pA u} 2, 12), 



(25) 



where f(u>) = l/[exp(o;/T) — 1] is the Bose-Einstein distribution function, and aj m = 2TrmT are the 
Matsubara frequencies. To extract the imaginary part, one uses the Dirac identity, Iml/(a; + ie) — 
—tt5(x), and analytic continuation, — iu) m —>■ uj + ie, for the retarded self-energy, 



ImS(w, k) 



N 1 



duii dui 2 [l + f(u>i) + f(u>2}] 8{uj — u>x — u> 2 ) 



/4Acr 

\~ J 2(2tt) 3 k J_ c 

dqi qi dq 2 q 2 6 (|c/i - q 2 \ < k < q x + q 2 ) p n (ui, qi)p n (u2, Hi) 



(26) 
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As discussed in Sec. [H] in the improved Hartree approximation of the O(N) model the spectral 
density of the pion is just a delta function (|16fl . 

(4Acr\ 2 A~ 1 f°° 
"ITT ) 7TT I ^1 du) 2 [l + /(^l) + /(wa)] S(u> - Ui - w 2 ) 

<2gi gi dg 2 92 © (|gi - 92! < k < gi + 92) [^(gi^fe)] -1 

x {<5[lji - Wtt^i)] 5[uj 2 - Un(qa)] + + u*(Qi)] ^2 + 

- S[ui - w^(gi)] <5[w 2 + w OT (g 2 )] - + w^gi)] 5[w 2 - w ff (g 2 )]} > (27) 



where w w (g) = \J q 2 + M% is the quasiparticle energy of the pion. To simplify this expression one 
uses the 5{lo ~ u>\ — ui 2 ) function to carry out the w 2 -integration, 



/4Act\ 2 N 1 r 

imE( ^ fc) = (— ) mil 



du)x[l + f(ui) + J{lo - uj\) 



poo 

x / dq 1 qi dq 2 g 2 6 (\qi - g 2 | < k < q 1 + q 2 ) [u 7! (q 1 )uj 7! (q 2 )Y 1 
Jo 

x {S[(a}! - w w (gi)] 5[w - wi - w,r(g 2 )] + + ^(gi)] 6[w -wi + w w (g 2 )] 

- - w^(gi)] 5[w - wi + w 7r (g 2 )] - + w OT (gi)] 5[w - a>i - u- K (q 2 )]} , (28) 

and the wi-integration with the help of the S[lu\ ± ^(gi)] -functions 

(4A(t\ 2 A^ 1 f°° 
—J J d QiQi d Q2q2 0(\q 1 -q 2 \<k<q 1 +q 2 )[u n (q 1 )Lu v (q 2 )]~ 1 

X ( {1 + f[(jJn{qi)] + f[oJ - 0J n (qi)]} 8[w - ^(gi) - u w (q 2 )] 

+{1 + f[-Vn(qi)] + /[w + w OT (gi)]} 5[cj + w w (gi) + c^vfe)] 
-{1 + /[w w (gi)] + /[w - w OT (gr)]} <5[w - w ff (gi) + w 7r (g 2 )] 

-{1 + /[-ww(?l)] + /[W + 5[W + UJ^(qi) - U^fa)]) ■ (29) 

The 5[cj + W7r(gi) + W7r(g 2 )] function has no support, because u> > and oj^{q) > 0. Note that evalu- 
ating the S[oj — u> v (qi) +u>^(q 2 )] and the 5[u> +u> w (qx) — w w (g 2 )] function, the latter two terms cancels 
each other. The remaining delta function can be transformed to a delta function in momentum- 
space, 5[w-w jr (gi)-a; w (g 2 )] = |w 7r (g )/go| [<%l - go) + <%i + go)] where g = \/P - ^-fe)] 2 - M 2 , 
is the root of the argument of the delta function. Note that in Eq. (|29(1 qi > and therefore there is 
no support of the 5(qi +qo) function. Carrying out the the gi-integration, using w — w w (go) = w w (gi), 
and relabelling q = q 2l one gets 

(4A<7 A 2 A^ 1 f 00 
— J g-^ / dq g[w 7r (g)] _1 6 (|g - g| < fc < g + g) 

x {1 + /K(g )]+/K(g)]}. (30) 

To calculate the limit fc — * 0, best one starts with Eq. I|29|l. uses the transformation 1 + /[w„-(gi)] + 
f[u - w*(?i)] - [1 - exp(-w/T)]{l + /K(gi)]}{l + /[w - w*(gi)]}, the identity 

,. © (|gi - 92| < k < gi +g 2 ) , , ql s 

hm =2 5(gi-g 2 ) (31) 

(to perform the gi-integration), and relabels q = q 2 , to get 

ImS(c,0) = ^y£[l-cxp(-^)]^ °dgg 2 %-2 W7r (g)][^(g)]- 1 

x {l + /[^(g)]}{l + /[w- W7r (g)]}. (32) 
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As discussed in detail in appendix E of Ref. [33j, this expression can be calculated analytically, 
using the Spence integral, 



Note that in |3^| the calculation is performed for the special case u> = m a , but this can be generalised 
without further problems. 

To calculate the real part of Eg. (|25fl . one has to integrate over the principal value of \/{loi+lo 2 —lo) 1 
denoted by P J . . . , 



/ &\rr\ 9 Af 1 /*°° 

[Re£(^)] 2 = (-) P f M 



1 + f(uj 1 ) + f(uj 2 ) 

2 ; 



poo 

X / dq 1 q 1 dq 2 q 2 Q (\qi - q 2 \ < k < qx + q 2 ) p 7T (oJ 1 ,q 1 )p 7T (uj2,q2)- (34) 
JO 

Again, one first uses the fact that the spectral density of the pion is just a delta function, which can 
be used to perform the uj\- and the w 2 -integration. Using trivial renormalisation, i.e., neglecting the 
divergent parts leads to 

(4A(t \ N 1 f°° 
~N~) 8tt2 k P J dqiqi dq2q2 e ~ 92 ' - k ~ qi +92 ' ) t^^ 1 )^^ 2 )] -1 

7[Wtt((?i)] + f[uJir{q 2 )} f[w*(qi)] +f[Uir{<l2)] 



oJ v (qi) + u> n (q 2 ) - oj uj n (qi) + Unfa) + u 

-/K(gi)] + /K(g 2 )] | -/K(gi)] + /K( g2 )] 
u}^(qi) - u w (q 2 ) - u> uj^qx) - oj w (q 2 ) + co 



(35) 



To evaluate the principal value numerically, in an appropriate way, one performs the following steps. 
First one introduces a new variable x = lo^qx) and transforms the q\ -integration to an rr-integration, 
with dq\ = u) v (qi) I qidx, and q\ — \J x 1 — M 2 = q* . Relabelling q = q 2 leads to 

/ A \ \ 2 AT 1 t°° 

[ReE(w,fe)] a = (-jfj J dqq[^{q)]- l dxQ (\q* - q\ < k < q* + q) (36) 

x \ f{x) + f[^(q)\ | /(.T) + /K(g)] | -I(x) + f[^{q)\ | -f(x) + f[u; 7T (q)} \ 
[ x + lo-x (q) — to x + uj 7! (q)+ui x — (^^(q) — U) x — (^^(q) + uj J 

Second one transforms the x-integration to a sum over xi , and uses the mean value theorem to put 
the terms with the distribution functions in front of the integrals 

[ReS( W ,fc)] 2 = \—j ^k P J dqqlco^qT 1 ^® (\q*~q\<k<q* + q) 

x ({f(Si) + /M?)]} dx ' 



x + U!„(q) - to 



-{f(x) + f[u,„(q)]} dx ' 



X + LUn (q) + UO 



+{-m + / 1+1 dx 1 



+ { - /t -'- , + / ^'^ l ^/ r ^ 1<fa; 3; -^ 1 ( g ) + a, J' 
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where Xi < x < Xi+i, and q* = yx — M 2 . Finally, the x-integrals can be performed analytically, 

—J ~./ d^K^r'E 9 (\q*-q\<k<q*+q) 

[u v (q) + x l+1 } 2 ~uj 2 



x {/(x) + /K(<z)]}ln 



+ {-/(£) + In 



w„.(g) + Xi] 2 - w 2 
w ff (g) - Xj+i] 2 - w 2 



[w ff (g) - x^ 2 - w 2 

In the limit k — > 0, using Eq. 1)31(1 . the g-integration can be carried out, 



tReE(w,0)] s 



/4,VV A 



V iV 7 4tt 2 
x ({/(£) + /K(g*)]} In 



i 



+ {-/(x) + /K(g*)]}ln 



^(g*) + x % ] 2 - lu 2 
^(g*) - x l+ i] 2 - 



[to w (q*) - Xi] 2 - uj 2 



(38) 



(39) 
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